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r <D - ABSTRACT 

: 

We derive constraints on the mass-temperature relation of galaxy clusters 
! from their observed luminosity-temperature relation and X-ray temperature func- 

tion. Adopting the isothermal gas in hydrostatic equilibrium embedded in the 
universal density profile of dark matter halos, we compute the X-ray luminos- 
ity for clusters as a function of their hosting halo mass. We find that in order 
to reproduce the two observational statistics, the mass-temperature relation is 

O fairly well constrained as T gas = (1.5 ~ 2.0) keV(M vir /10 14 /if 1 M o ) a5 ~ - 55 , and a 

simple self-similar evolution model (T gas oc M^ 3 ) is strongly disfavored. In the 
cosmological model that we assume (a ACDM universe with Q = 0.3, Ao = 0.7 
and h 7Q = 1), the derived mass-temperature relation suggests that the mass 

"Sa '. fluctuation amplitude a% is 0.7-0.8. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general 
^ \ — X-rays: galaxies 

d; 

1. Introduction 

While clusters of galaxies are relatively simple dynamical systems that consist of dark 
matter, stars, and X-ray-emitting hot gas, their thermal evolution is not yet fully under- 
stood. This is clearly illustrated by the well-known inconsistency of the observed X-ray 
luminosity-temperature (Lx-T) relation, Lx oc T 3 (e.g., David et al. 1993; Markevitch 1998; 
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Arnaud & Evrard 1999) against the simple self-similar prediction L x oc T 2 (Kaiser 1986). 
Conventionally this is interpreted as evidence for preheating of intracluster gas; additional 
heating tends to increase the temperature and the core size of the cluster, and to decrease 
the central density and the luminosity. Since the effect is stronger for less massive systems, 
the slope of Lx-T relation becomes steeper than that of the self-similar prediction (Evrard 
& Henry 1991; Kaiser 1991). 

In addition, the mass-temperature (M-T) relation of clusters is also poorly determined. 
Although it is conventionally assumed that the gas shock heating is efficient enough and 
the temperature of the intracluster gas reaches the corresponding virial temperature of the 
hosting halos, this should be regarded as a simple working hypothesis. Nevertheless, the 
cosmological parameters derived from the cluster abundances are sensitive to the adopted 
M-T relation. While this has already been recognized for some time (e.g., Figs.5d and 6d 
of Kitayama & Suto 1997), Seljak (2002) recently showed in a quantitative manner that the 
use of the observed M-T relation by Finoguenov, Reiprich, & Bohringer (2001) decreases 
the value of the mass fluctuation amplitude at 8 h^ 1 Mpc, as, by ~ 20% where h is the 
Hubble constant Hq in units of 100 km s _1 Mpc -1 [we use, however, the dimensionless 
Hubble constant h-jQ = H /(70 km s _1 Mpc" 1 ) in the following analysis]. Therefore, the 
independent derivation of the cluster M-T relation is important both in understanding the 
thermal history of the intracluster gas and in determining the cosmological parameters. 

Our primary aim in this paper is to find the M-T relation of clusters that reproduces 
the observed Lx-T relation and X-ray temperature function (XTF). The reason we focus on 
M-T relation is as follows: since recent iV-body simulations strongly indicate the universality 
of the density profile of the hosting halos of clusters (Navarro, Frenk, & White 1996; Moore 
et al. 1998; Jing & Suto 2000), the intracluster gas density profile in hydrostatic equilibrium 
with the underlying dark matter can be computed (Makino, Sasaki, & Suto 1998; Suto, 
Sasaki, & Makino 1998) for a given mass of the halo. This enables one to make a reliable 
prediction for the Lx-T relation once the M-T relation is specified. In turn, one can obtain 
the M-T relation that reproduces the observed Lx-T relation without assuming an ad hoc 
model for the thermal evolution of intracluster gas. 

In what follows, we parameterize the M-T relation as a single power law, and derive the 
best-fit values of their amplitude and slope from the observed Lx-T relation and the XTF. 
The result is compared with the recent observational studies by Finoguenov et al. (2001) and 
Allen, Schmidt, & Fabian (2001). We also discuss the implications for the value of cr 8 from 
cluster abundances. Throughout the paper, we adopt a conventional ACDM model with 
density parameter f2 = 0.3, cosmological constant A = 0.7, dimensionless Hubble constant 
h 70 = 1, and baryon density parameter Qb = 0.04 h^Q. 
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2. Model of intracluster gas in dark matter halos 

We first outline the model of the dark matter halo and the isothermal gas density profile 
embedded in the halo, which are essential in predicting the X-ray luminosity of clusters as a 
function of the mass of the hosting halo. 



2.1. Dark Matter Density Profile 

We adopt the specific density profile of dark matter halos of mass M vir , given as 

p(z) 5 c (M vir ) 



r < r vir , 



Phaio(r;M vir ) = { (r/r s r(l + r/r s y-« ™' (1) 

r > r vir , 

where p(z) = fioPeo(l + z) 3 is the mean density of the universe at z, p c o is the present critical 
density, <5 c (M vir ) is the characteristic density excess, and r vir and r s are the virial radius and 
the scale radius of the halo, respectively. In practice, we focus on two specific profiles: a = 1 
(Navarro et al. 1996) and a = 3/2, indicated by higher resolution simulations (Moore et al. 
1998; Jing k Suto 2000; Fukushige & Makino 2001). 



(M^(^F-) ■ (2) 



The virial radius r vir is defined according to the spherical collapse model as 

1/3 

V47rpA n i 

and the approximation for the critical overdensity A n i = A n i(fi ,Ao) can De found in Ki- 
tayama & Suto (1996). The two parameters r s and r vir are related via the concentration 
parameter, 

c = ciM '- z) = TjM-W (3) 

In the case of a — 1, we use an approximate fitting function with the same functional 
form as that of Bullock et al. (2001), 

/ !\ it \ c norm ( -^vir \ / ,\ 

Cfl( M vl „ z) = — ^___ r _ j . (4) 

Since the c-M vir relation has a fairly large intrinsic scatter, we adopt the same value of the 
power-law index (—0.13) as indicated by Bullock et al. (2001), but perform our own fit to 
their Figure 4 to determine the value of the coefficient c norm . In fact, the original coefficient 
given by Bullock et al. (2001) does not seem to fit their data, and we find c norm = 8^7- 
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The quoted errors do not represent the fitting error to the mean relation, but correspond to 
±1(7 for the intrinsic distribution around the mean relation. The uncertainty with respect 
to the adopted power-law index (—0.13) is effectively included in the above quoted errors 
for c norm . For a ^ 1, we rescale the amplitude of the concentration parameter according to 
Keeton & Madau (2001) as c(M vir ,z) = (2 - a)c B (M vil , z). 

The condition that the total mass inside r vir is equal to M vir relates 5 C to c as 



5 c (M vil ) 



^nl C 



3 m(c)' 



(5) 



where 



m{x) 



ln(l + x) -x/(l + x) 



ln( dx + y/l+x) - 



x 



1 + x 



a = 1, 
a = 3/2. 



(6) 



2.2. Hydrostatic Equilibrium Gas Distribution 

We further assume that the intracluster gas is isothermal and in hydrostatic equilibrium, 
which is a reasonable physical approximation. Under the gravitational potential of the 
above dark matter halos, the isothermal gas density profiles in hydrostatic equilibrium are 
computed analytically as (Suto et al. 1998) 



where 



and 



Pgas(r) = p g as,o exp[-5/(r/r s )], (7) 

B = m{c) T gas (8) 
1 — — ln(l + ar) a = 1, 

/(*) = < [TTx- 2 ^ 

2d ln(v^+ y/l + x) a = 3/2. 

V Ob Ob 

In equation (8), T vir is the virial temperature, which we define as 

A: B T vir = \pmf^ oc M% 3 , (10) 

where k-Q is the Boltzmann constant, G is the gravitational constant, and fim p is the mean 
molecular weight. We adopt /i = 0.6 assuming that the gas is almost fully ionized with the 
mass fractions of helium Y = 0.24 and metals Z = 0.3 Z & (Zq = 0.02). On the other hand, 
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the gas temperature T gas is determined from the M-T relation, as described in detail in §2.3. 
The central gas density p gaSj o is computed so that 

p gas (r)A7rr 2 dr = / gas M vir (^) , (11) 

where / gas is the hot gas fraction of baryon mass in the cluster described in §2.4. The X-ray 
luminosity of clusters is computed as 

rr v i T 

L x = 4vr / A(T gas , Z) 
Jo 

In practice, we adopt the bolometric cooling function of Sutherland & Dopita (1993) for 
A(T gas ,Z). 



z 



^gasV 



r dr. 



(12) 



2.3. Mass-Temperature Relation 

As emphasized in the above, we are mainly interested in the M-T relation of clusters. 
Most of previous theoretical studies adopted the self-similar relation for the M-T relation: 

2 /3 

Tgas — T v - ir oc M vir . Recent observations, however, indicate a departure from this relation. 
Finoguenov et al. (2001), for instance, obtained 

/ M \ 0.54±0.02 

r » = (M3 ± °' 07) keV {w^M-J • (13) 

where Ma c is the total mass enclosed within the radius r& c at which the mean interior density 
is A c times the critical density of the universe, and T cw is the emission-weighted temperature. 
Finoguenov et al. (2001) estimated M 50 o from the observed X-ray luminosity density profile 
assuming that the gas is in hydrostatic equilibrium. Similarly Allen et al. (2001) found 

(yr \ 0.65±0.09 

— — — — j— ] , (14) 

where T 250 o is the gas-mass-weighted temperature within r 2 5oo (see also Ettori, De Grandi, 
& Molendi 2002). 

In order to generalize the various choices in the M-T relation, we adopt the parameter- 
ization 

M ■ x PMT 



T ^ M ^ =T ^[w^wJ ■ (15) 

For reference, the self-similar model T gas = T vir with equation (10) corresponds to (T gaS)0 ,pMT) 
(1.1 keV, |). Since M 500 and M 250 o quoted in equations (13) and (14) are different from M vir , 
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we translate them to M vir , properly taking into account the dark matter density profiles as 
described in Appendix A (the relations between M& c and M vir for several different values 
for the overdensity A c are plotted in Fig. 9 below). In doing so, we also convert the cluster 
data at the individual redshifts into those at z — 0, taking into account the difference of 
the Hubble parameter at z and Hq for our assumed cosmology (Qq = 0.3, Ao = 0.7, and 
h 70 = 1.0). Then we find that the observed samples are well fitted to the relations 



0.54±0.02 



(1.92 ± 0.06) keV ( — ) a = l, 

r g .,H y° ll $ M °L^ tie) 
(L88±0 - 06)keV (ib^f^) a ^ 3A 



for Finoguenov et al. (2001), and 



M vir 



0.57±0.12 



'1.53 ± 0.56) keV - T a = l, 

;L45±0 - 54)keV (l^|l5e) °^ /2 ' 



for Allen et al. (2001). 



2.4. Hot Gas Mass Fraction 

Hot gas fraction / gas also plays a central role in predicting the X-ray luminosity of 
clusters. In many theoretical analyses, it is often assumed that / gas is independent of the mass 
of the hosting halos just for simplicity. Of course, this is not a good approximation because 
the gas in less massive systems is expected to have cooled more efficiently at high redshifts 
on average, and thus / gas should be a monotonically increasing function of the halo mass. 
This qualitative feature is supported by both observations and hydrodynamical simulations. 
Mohr, Mathiesen, & Evrard (1999), for instance, measured the gas mass fraction within 
as a function of gas temperature in the range 3 keV < T gas < 10 keV. Extrapolating their 
result to lower temperatures, and assuming for simplicity that the gas mass fraction at the 
virial radius is equal to that at r 50Q , we obtain 

/ gas = min 

This is our fiducial model for the hot gas fraction in the present analysis, and for comparison, 
we also consider a simple model / gas = 0.8, in which the gas mass fraction is independent of 
halo mass and gas temperature. 



0.92/i 



-3/2 
70 



T 

-L fy. 



6 keV 



0.34 



(18) 
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Strictly speaking, we have to convert the value of the gas fraction of Mohr et al. (1999) 
defined at r 50Q to that at the virial radius. In practice, however, the observational error of the 
value is significantly larger than the difference of the conversion, and thus we assume here 
that / gas (^5oo) — igas( r vir)- We put an additional condition that the hot gas to baryon fraction 
in clusters inside their virial radius does not exceed the cosmological baryon fraction for our 
assumed cosmological parameters (Oo = 0.3 and Qb = 0.04 /if 2 ). In fact, the observed gas 
to dark matter fraction has a large scatter (see Fig. 14 of Mohr et al. 1999) and is consistent 
with the upper bound that we set here. Actually, if we adopt the observational M-T relation 
of Finoguenov et al. (2001), we can constrain the gas mass fraction so as to reproduce the 
observed L x -T relation and XTF. As shown in Appendix B, the resulting constraint is fairly 
consistent with equation (18). 



3. Mass-temperature relation 

3.1. Constraints from the Luminosity- Temperature Relation 

As mentioned above, the simple self-similar model prediction L x oc T 2 is too shallow to 
be consistent with the observation (Lx oc T 3 ). This means that heating/cooling processes 
in addition to the shock heating are important in the thermal evolution of intracluster gas. 
Apart from the physical mechanism of the additional thermal processes, there are three 
possibilities that might modify the mass dependence of X-ray luminosity (see eq. [12]) and 
steepen the resulting L x -T relation. First, the gas density profile may be significantly flatter 
for less massive systems. Second, the mass dependence of the hot gas mass fraction is strong 
as /g as oc M^/ r 3 . Finally, the mass-temperature relation is T gas oc M^/ 5 . In practice, a realistic 
model should be a combination of those three effects to some extent. The gas density profile 
can be specified completely from the hydrostatic equilibrium assumption. As for the hot gas 
mass fraction, we adopt the observed relation as well as the simple constant fraction. Then 
the model described in the previous section enables one to compute the X-ray luminosity of 
a given halo, and to examine whether the observed L x -T relation can be reproduced from 
observed M-T relation and the gas mass fraction. 

Before doing so, let us first look at the L x -T relation predicted from the observed M-T 
relations (eqs. [16] and [17]) and gas mass fraction (eq. [18]) combined with the isothermal 
gas density profile (eq. [7]). Figure 1 compares those predictions against 52 X-ray clusters 
with temperature higher than 2.5 keV from the sample of Ikebe et al. (2002) (excluding the 
two clusters with no reliable estimates for their temperatures) . For comparison, we also plot 
the result for the self-similar model (eq. [10]). Incidentally, we performed all the analysis 
below both for a = 1 and 3/2, but their difference turns out to be very small. Thus, we 
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show the results for a = 1 except for the combined contour plot (see Fig. 8 below). 

This clearly illustrates that the predicted L x -T relation is very sensitive to the assumed 
M-T relation; as is well known, the simple self-similar model (Fig. 1, dashed line) is incon- 
sistent with the observations by a wide margin. If the halo density profile is well described 
by equation (1) with 1 < a < 3/2, the M-T relation of Allen et al. (2001; Fig. 1, dotted line) 
is in good agreement and that of Finoguenov et al. (2001; solid line) leads to an acceptable 
result. This conclusion in turn indicates that the Lx-T relation provides a good diagnosis of 
the underlying M-T relation, which is as yet poorly determined observationally. 

Therefore, we next attempt to find the range of parameters (pmt and T gaSj o in eq. [15]) 
which reproduces the observed Lx-T relation. As is clear from Figure 1, the observed data 
have intrinsic dispersions. Thus, we first divide the cluster sample in nine temperature bins 
so that each bin contains five or six clusters. Then we compute the mean temperature, and 
the mean luminosity and the standard deviation of clusters in each bin. The results are 
plotted by open circles with error bars for the luminosity in Figure 2. Then we perform a 
X 2 fit to the binned data. 

The result is plotted as confidence contours on the T gaSi o-pMT plane in Figure 3. We 
estimate the relative confidence levels with respect to the best-fit values assuming that 

AX 2 = X 2 (T gas , ,PMT) ~ X 2 (^gas,0,min,PMT,min) (19) 

follows the x 2 distribution of the 2 degrees of freedom (e.g., Press et al. 1992, chap. 15), 
where T gaSi o, m m and pmt, mm are the best-fit parameters that minimize the value of x 2 ■ Upper 
and lower panels correspond to the gas mass fraction of / gas = 0.8 and of equation (18) from 
Mohr et al. (1999). The three contour curves represent the 1, 2, and 3 a confidence levels 
derived from the value of Ax 2 - Crosses show the positions of the best-fit values indicated 
in each panel. For comparison, we plot the self-similar model prediction, and observational 
estimates by Allen et al. (2001) and Finoguenov et al. (2001) as open triangles, open circles, 
and filled circles, respectively. Figure 3 rephrases the visual impression from Figure 1 in 
a more quantitative way; the M-T relation of Allen et al. (2001) is inside the 1 a contour 
and that of Finoguenov et al. (2001) is located just outside the 3 a confidence level, and 
marginally consistent within the large error-bars. While the concentration parameter of 
dark matter halos has a fairly broad distribution corresponding to c norm = 8^27 * n ec L ua tion 
(4), it does not lead to any significant difference (compare dotted and dot-dashed lines with 
solid lines in Fig. 3). Thus, we fix the proportional constant of the concentration parameter 
Cnorm = 8 in what follows. 

The resulting Lx-T and M-T relations for the best-fit parameters are plotted in Figures 2 
and 4. While our best-fit models actually reproduce the observed Lx-T relation (the best-fit 
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value of x 2 P er degree of freedom is shown in each panel of Fig. 3), they also seem to be in 
reasonable agreement with the observed M-T relation. 



3.2. Constraints from X-ray Temperature Function 

One can also infer the empirical shape of the M-T relation from the requirement that 
it reproduces the observed XTF of clusters. Comparing with that from the L x -T relation, 
this methodology provides a fairly (even if not entirely) independent constraint on the M-T 
relation, mainly in two important aspects: (1) the model dependence on the density profile 
enters only through the flux limit, Sn m , of the observed sample, and (2) the prediction, on 
the other hand, is sensitive to the adopted mass function of dark matter halos, and therefore 
to the value of ag in particular (remember that we fix the other cosmological parameters in 
the present analysis). 

For this purpose, we again use the 54 X-ray clusters with temperature higher than 
2.5 keV from the sample of Ikebe et al. (2002). The corresponding flux limit is Su m = 2 x 10~ n 
erg s -1 cm -2 in the (0.1-2.4) keV band, and the total sky coverage is 8.14 sr. For definiteness, 
we adopt the gas mass fraction given by equation (18). As for the halo mass function, 
we adopt both an analytic model by Press & Schechter (1974) and a fitting model to the 
numerical simulations by Jenkins et al. (2001), 

dn p m d J n f\ (20) 



d\nM viT M vir v ' dh\M, __ 

f(a) = O.SlSexp^-lhKT^ + O.Gll 3 ' 8 ) . (21) 

We define the rms variance of linear density field a as 

a 2 (M) = 4tt / P(k)——^[sm(kR) - kRcos(kR)]k 2 dk, (22) 
J (kR) 

where R = (SM/Anp) 1 ^ 3 and P(k) is the linear power spectrum. We perform the x 2 fit to 
the data with respect to the three free parameters: pmt and T gaS)0 in equation (15), and a$. 

In order to avoid uncertainties in constructing the conventional XTF data (e.g., the 
definition of V max as discussed by Ikebe et al. 2002), we directly use the number count of 
clusters with the flux limit S\ iTa per unit solid angle of the sky, appropriated binned according 
to their temperatures. In practice, we divide both the data and our predictions into five 
temperature bins (2.5-3.6, 3.6-5.0, 5.0-6.4, 6.4-8.0, and >8.0 keV), so that the errors do 
not correlate with one another, and perform the x 2 fit- We assign the Poisson error to the 
number count in each bin. 
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The result is plotted as confidence contours on the T gaS)0 -pMT plane in Figure 5. Again 
we estimate the relative confidence levels with respect to the best-fit values assuming that 

^X =X (^gas,0) PMT, °"8,local min) — X (^gas,0,min, PMT,min, °8 ,min ) (23) 

follows the x 2 distribution of the 2 degrees of freedom, where T gaSi0mm , Pmt, min an d o"s,min 
are the best-fit parameters that minimize the value of x 2 , and og.iocai min is the value of a 8 
minimizing the x 2 with given T gaS) o, and Pmt- Upper and lower panels of Figure 5 correspond 
to the mass functions of Press & Schechter (1974) and Jenkins et al. (2001), respectively. 
The three solid contour curves represent the 1, 2, and 3 a confidence levels derived from 
the x 2 fi^ with the filled squares indicating the points of the highest significance. Dotted 
contours, on the other hand, indicate the best-fit value of <r 8 at the given location on the 
Tgasfi-puT plane. The self-similar model prediction, and observational estimates by Allen et 
al. (2001) and Finoguenov et al. (2001), are plotted by open triangles, open circles, and filled 
circles, respectively. 

The upper panel in Figure 6 illustrates that the best-fit model reproduces the observed 
XTF nicely. Again just for comparison, we plot the self-similar model predictions, which do 
not fit the observed XTF at all. Moreover, the lower panel indicates that the degree of fit is 
indeed sensitive to the value of cr 8 ; the best-fit value is 0.7-0.8. 

Although the above conclusion may seem inconsistent with the previous claims (Viana 
& Liddle 1996; Eke, Cole, & Frenk 1996; Kitayama & Suto 1997) that the self-similar model 
does fit the observed XTF with a 8 = 0.9-1.0, it can be explained by the following reasons. 
First, the previous XTF data had much larger errors and thus still allowed a wider range 
of theoretical models. Second, the latest XTF data that we adopted have a systematically 
smaller amplitude at T gas > 5 keV than the previous ones (see Fig. 6 of Ikebe et al. 2002). 
Third, the adopted Lx-T relation is different; our current self-similar model corresponds to 
Lx oc Tg as , while the previous analyses used (sometimes implicitly) Lx oc T| as to construct 
the "observed" XTF data via the conventional V max method (e.g., Eke et al. 1996). Finally, 
the recent mass function of Jenkins et al. (2001) predicts more massive halos than that of 
Press & Schechter (1974), which has been widely used in the previous analyses. 

Figure 7 demonstrates the above features. For a direct comparison with the analyses 
of Kitayama & Suto (1997), we adopt T gas = 0.8T vir and the L x -T relation from equation 
(3) of Kitayama & Suto (1997), with their fiducial choice for the other parameters. The 
goodness of the fit turns out to be significantly degraded compared to the previous result, 
mainly because of the systematically smaller amplitude of the XTF, as well as reduced error 
bars. Nevertheless, this methodology fully reproduces the best-fit value of erg = 0.9-1.0, as 
in the previous one. 
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In fact, the combination of those effects also explains why the reanalysis of the cluster 
abundance by Seljak (2002) yielded <j 8 ~ 0.7 in the standard ACDM model when adopting 
the M-T relation of Finoguenov et al. (2001). Incidentally, the smaller value of <7g seems 
consistent with the recent joint analysis of cosmic microwave background and large-scale 
structure (Efstathiou et al. 2002), which suggests as = 0.6-0.7. 

4. Discussion and conclusions 

We have presented the constraints on the (empirically parameterized) mass-temperature 
relation of galaxy clusters from the two different observational data, the luminosity-temperature 
relation and the temperature function. 

We summarize in Figure 8 our constraints on the T gaSi0 -pMT plane, adopting the observed 
gas mass fraction (eq. [18]). The fact that the simple self-similar evolution model T gas oc M^/ 3 
fails to explain the observed Lx-T relation is well known and not at all new. Rather, it should 
be noted that the mass-temperature relation of galaxy clusters is fairly well constrained 
by combining the observed Lx-T relation and XTF. It is encouraging that T gas = (1.5 ~ 
2.0) keV(M vir /10 14 M ) a5 ~°- 55 barely satisfy the two constraints simultaneously. This 
conclusion applies for 1 < a < 3/2 as long as the dark halo density profile is described by 
equation (1). 

In addition, our analysis implies that the amplitude of the mass variance <r 8 in the 
standard ACDM model should be 0.7-0.8. These values are significantly smaller than the 
previous estimates (Viana & Liddle 1996; Eke et al. 1996; Kitayama & Suto 1997) but in 
better agreement with more recent results (Seljak 2002; Efstathiou et al. 2002). 

Since our current analysis has adopted a simple parameterized model for the mass- 
temperature relation, the result should be understood by a physical model of the thermal 
evolution of the intracluster gas. For instance, our result may be qualitatively explained by 
a kind of phenomenological heating of T gas (M vir ) = T vir (M vir ) + 1 keV. We plan to examine 
the implications for the possible heating sources from the derived mass-temperature relation 
of galaxy clusters using the Monte-Carlo modeling of merger trees. 

We thank an anonymous referee and Alexis Finoguenov for useful comments. This 
research was supported in part by the Grant-in-Aid for Scientific Research of JSPS (12640231, 
14102004). S.S. gratefully acknowledges support from TMU President's Research Fund for 
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A. Mass-radius relation at a given overdensity 

In order to determine the mass of a cluster, one has to specify the radius from its center. 
Observationally, this is often set by the specific fractional overdensity A c with respect to the 
critical density of the universe, p c . Thus the resulting mass M& c is written in terms of the 
corresponding radius r& c as 

M Ac = ^7rri c A cPc . (Al) 

Most theoretical studies, on the contrary, usually define the mass of dark matter halo M v ; r at 
its virial radius r vir , adopting the spherical nonlinear collapse model, i.e., A c = fi A n i(fi , A ) 
(see Kitayama & Suto 1996 for details). 

Once the density profile of dark matter halos is specified, Ma c can be easily translated 
into M vir . For the profile that we adopted in this paper (eq. [1]), mass inside a radius r is 
written as 

M(r) = 47r5 c p(z)r 3 s m(r/r s ), (A2) 

where m(x) is defined in equation (6). Thus, if A c is given, one can compute r^ c and 
M/s. c /M vir from the following relation 

M Ac = m(r A Jr s ) = r^A^p, 
M vir m(c) r vir A n i p' 

Figure 9 plots MA C /M vir {upper panels) and r Ar (lower panels) for several choices of A c . 



B. Hot gas mass fraction 

The strategy that we have adopted in the present paper is to determine the permitted 
parameter region of the M-T relation of clusters from the L x -T relation and XTF, assum- 
ing a specific model for the hot gas mass fraction f gas (M y i T ) inspired by the observation. 
In principle, however, we may repeat the similar procedure and derive the constraints on 
/gas (M v i r ), assuming the observed M-T relation instead. Since the observational uncertainty 
for the hot gas mass fraction (eq. [18]) that we adopted is fairly large, this approach is useful 
in understanding the dependence of our conclusion on the model for the gas mass fraction. 

For this purpose, we also parameterize the gas mass fraction by a single power law, 

(rp \ Pgas 

TkeV J > (Bl) 

and perform the x 2 fit to the L x -T relation varying the two free parameters, / gas ,o an d p gas - 
We repeat the similar fitting XTF with / ga s,o, Pgas and erg. 
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The result is plotted in Figure 10, where we use the observed M-T relation of Finoguenov 
et al. (2001). Solid contour curves represent the 1, 2, and 3 a confidence levels derived from 
the x 2 fit to the Lx-T relation, while dotted curves show the 1, 2, and 3 a confidence 
levels from the XTF. Clearly, both constraints again are simultaneously satisfied for / gas = 
(0.5 ± 0.1)(T gas /l keV)°' 4±015 , and in fact agree with the observed gas mass fraction of Mohr 
et al. (1999). We also plot the condition that / gas (10 keV) = 1 by a dashed line, so that the 
hot gas fraction does not exceed the cosmological average of the baryon fraction (assuming 
that £l B = 0.04 hjQ and f2 = 0.3). If this condition should be satisfied, the acceptable 
parameter range almost exactly corresponds to the observed gas mass fraction of Mohr et 
al. (1999). 



-16- 



10 2 r 



\ 
SO 

u 

cv? 

o 



10 - 



1 - 



10 



-1 



f : Mohr et al. w — - 

gas ^ 




T . 

vi r 

Allen et al. 

Finoguenov et al. 



4 



gas 



6 

[keV] 



8 10 



Fig. 1. — X-ray Lx-T relations derived from the observed M-T relations (solid line, eq. [16]; 
dotted line, eq. [17]) assuming the observed gas mass fraction (eq. [18]). A simple self-similar 
model prediction with T gas = T vir is plotted by the dashed line for comparison. The data 
with error bars indicate 52 X-ray clusters with temperature higher than 2.5 keV from the 
sample of Ikebe et al. (2002). 
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Fig. 2. — Best-fit Lx-T relations from the parameterized M-T relations (crosses in Fig. 3). 
The nine open circles with error bars indicate the binned data described in the text. Solid 
and dotted lines correspond to the cases assuming the gas mass fraction of / gas = 0.8 and of 
equation (18), respectively. 
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Fig. 3. — Constraints on the parameterized M-T relation of clusters from the y 2 fit to the 
binned Lx-T data. Upper and lower panels adopt the gas mass fraction of / gas = 0.8 and 
of eq. (18), respectively. The three contour curves represent the 1 cr(68.4%), 2 cr(95.6%), 
and 3 cr(99.7%) confidence levels for c norm = 8. Crosses show the positions of the best- 
fit values indicated in each panel. For comparison, also plotted are the self-similar model 
prediction (open triangles), and observational estimates by Allen et al. (2001; open circles) 
and Finoguenov et al. (2001; filled circles). Dotted and dot-dashed contours indicate the 
results of 3 <r(99.7%) confidence levels for c nor m = 10 and 5.3, corresponding to ±1 o deviation 
from the average. 
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Fig. 4. — Predicted M-T relations compared with observations. Solid and dotted lines 
indicate the M-T relations derived from the observed Lx-T relation for / gas = 0.8 and 
eq. (18), respectively (corresponding to the crosses in Fig. 3). The observational data points 
and the fits are taken from Finoguenov et al. (2001; dots with error bars) and Allen et 
al. (2001; six filled circles with error bars) after correction for the difference of the mass 
definitions (cf. Appendix A). A simple self-similar model prediction with T gas = T vir is 
plotted in dot-dashed line for comparison. 
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Fig. 5. — Constraints on the parameterized M-T relation of clusters from the \ 2 fit to the 
observed XTF. Upper and lower panels adopt the mass functions of Press & Schechter (1974) 
and Jenkins et al. (2001), respectively. The three solid contour curves represent the 1, 2, and 
3 a confidence levels, with the filled squares indicating the points of the highest significance. 
Dotted contour curves show the best-fit values for <j$ (which label each curve). Symbols are 
as in Fig. 3. 
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Fig. 6. — Best-fit XTF (upper panel) from the parameterized M-T relation compared to the 
cluster sample of Ikebe et al. (2002) and the corresponding \ 2 as a function of a 8 (lower 
panel). Solid and dotted curves adopt the mass functions of Jenkins et al. (2001) and Press 
& Schechter (1974), respectively. A simple self-similar model prediction with T gas = T vir is 
also plotted just for illustration. 



-22 - 



10 =- 



GO 



CO 

A 



us 



E- 

A 



10" 1 =- 



n 1 — i — i — i — i — i — | — 

observed LT 




W keV ]= - 88 

p MT =2/3 



J I I I I L 




T gas [keV] 



20 







- 1 \ 1 ' 

T =0.8T . \ '. / 

gas vir \ , / 


1 1 1 1 

A X z=4.0 - 




Ax 2 =1.0 - 






Ax 2 =4.0 " 




Ax 2 =1.0 " 


Jenkins 

PS 


best fit a e 
0.91 (J) 
1.00 (PS) " 


observed LT 

i 1 , 1 , 


I.I." 



0.7 0.8 



0.9 



1.1 



1.2 



C7 C 



Fig. 7. — Same as Fig. 6, but for a self-similar model adopting the observed Ly^-T relation 
and a slightly smaller proportional factor (T gas = 0.8T vir ). 
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Fig. 8. — Joint constraints on the parameterized M-T relation from the binned Lx-T data 
(solid curves; Fig. 3) and from the XTF (dotted curves; Fig. 5). Upper and lower panels 
adopt the mass functions of Press & Schechter (1974) and Jenkins et al. (2001), respectively, 
a = 1 (Left) and a = 3/2 (Right). Symbols are as in Figs. 3 and 5. 
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Fig. 9. — Ratio, M^ c /M vir (upper panels) and radius, r& c (lower panels) , for a given overden- 
sity A c as a function of M vir , a = 1 (left) and a = 3/2 (right). Solid, dotted, short-dashed, 
and long-dashed lines correspond to the overdensities A c = Q A n \, 200, 500, and 2500, re- 
spectively. The standard critical overdensity predicted in the nonlinear spherical collapse 
model Qq A n \ is ~ 100 in our fiducial values of the cosmological parameters with Qq — 0.3 
and Aq = 0.7. 
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Fig. 10. — Joint constraints on the parameterized hot gas mass fraction relation from the 
observed Lx-T relation (solid curves) and from the XTF (dotted curves). The fit to the 
observational data by Mohr et al. (1999) is indicated as an open circle with error bars. 
Dashed lines indicate the condition that / gas (10 keV) = 1. 



